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I. INTRODUCTION TO CLUSTER-CLUSTER 
AGGREGATION (CCA) 

Cluster-cluster aggregation (CCA) refers to the statis- 
tical dynamics of a large collection of particles (or clus- 
ters) which move around under the action of some trans- 
port process and which, in addition, have a probability 
to coalesce irreversibly when they collide. The transport 
is usually stochastic, diffusion being a commonly stud- 
ied case, and coalescence of clusters conserves mass. In 
general, both the rate of coalescence and the transport 
coefficient are functions of cluster mass. Such dynamics 
are relevant to many physical phenomena. Commonly 
cited applications include the study of aerosols [T], the 
development of structure in the expanding universe [2], 
clustering of algae [31 and cloud droplet coalescence [3]. 
From a theoretical perspective, CCA can be modelled us- 
ing relatively simple interacting particle systems and can 
thus provide a testing ground for statistical mechanics 
theories of complex systems. The theoretical objective 
is to understand the statistical properties of the mass 
distribution of the clusters. 

CCA is an intrinsically non-equilibrium process. Ag- 
gregation is irreversible so there is no possibility for the 
statistical dynamics to have detailed balance. In some 
circumstances it may be appropriate to consider frag- 
mentation of clusters in addition to aggregation. In such 
cases, detailed balance may be achievable. We do not 
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consider such situations here however. Two different sce- 
narios may be studied which we refer to as the decay and 
forced cases. In the decay case, the more common of the 
two, the system is initially populated with a large num- 
ber of small particles (monomers) which are considered to 
be the fundamental unit for the aggregation process and 
whose mass can accordingly be taken to be one. The size 
distribution subsequently spreads in the space of cluster 
masses as larger clusters are generated by the aggrega- 
tion of smaller ones. This time evolution proceeds for 
all time (we do not consider models exhibiting a gelation 
transition here) often in a self-similar fashion. In the 
forced case, the system initially has no clusters but we 
continually deposit monomers into the system at a con- 
stant rate from an external source. The size distribution 
again spreads in the mass space but, due to the constant 
replenishment of the monomers by the external source, 
eventually reaches a stationary state at any given cluster 
size. In this stationary state, the depletion of clusters 
of a given size by aggregation events is exactly balanced 
by the formation of such clusters by the aggregation of 
smaller ones. This results in a constant flux of mass in 
the mass space. An analogy may be drawn with the flux 
of energy through scales in the Kolmogorov cascade of 
hydrodynamic turbulence '5 . We shall be concerned pri- 
marily with the stationary forced case. 

The most basic statistical measure is the average mass 
density, N{m, t), giving the average density of clusters of 
mass m at time t. At the mean-field level, N{m, t) sat- 
isfies the Smoluchowski equation [6]. For many applica- 
tions, the collision kernel in the Smoluchowski equation, 
which reflects the mass dependence of the underlying mi- 
croscopic coalescence rates, is a homogeneous function 
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whose degree we shall denote by For a broad class of 
such kernels, the Smoluchowski equation has self-similar 
solutions which describe scaling behaviour of N{m,t) in 
both m and t. In particular, N{m,t) has power law 
asymptotics for N(m, t) as a function of mass. In the de- 
cay case, we often find N{m, t) ^ m^*^^+^). In the forced 
case, the stationary state is given by N{m) ~ 
[51 [7]. The reviews by by Aldous and Leyvraz [5] 
provide comprehensive discussions of the Smoluchowski 
equation and its scaling solutions. 

One of the most interesting aspects of CCA is that this 
mean field description can fail under certain conditions. 
This results in a transition to a fluctuation dominated 
regime with completely different scaling properties for 
the large m asymptotics. The fundamental assumption 
underpinning mean field theory is the absence of spa- 
tial correlations. In particular, the joint distribution of 
clusters having masses mi and m2 at a given point, x, 
which is linked to the probability of collision, factorises 
as the product of the one-point distributions of clusters 
with masses mi and TO2- Likewise for higher order joint 
distributions. This works well in the so-called reaction 
limited regime when the density of clusters is sufficiently 
high that each aggregating cluster has a large number 
of potential partners to with which to aggregate. In 
the most commonly considered case, in which the trans- 
port process is simple diffusion, this assumption turns 
out to be inconsistent for large masses if the spatial di- 
mension is less than or equal to two. The reason is that 
in two dimensions and below random walks are recur- 
rent so that clusters which encounter each other once 
are highly probable to do so infinitely many times and 
hence to aggregate. Heavy clusters thereby generate an 
effective depletion zone around themselves. The dynam- 
ics thus becomes diffusion limited and dominated by dif- 
fusive fluctuations making the system much more diffi- 
cult to analyse. The only case for which a systematic 
study of the diffusion limited regime has been done is 
the case of constant kernel in which the aggregation rate 
is independent of mass (C = 0). For the forced case in 
d < 2, the diffusion limited stationary mass density can 
be shown to scale as N(m) ^ 77i-(2'i+2)/(<i+2) ^ifhich is 
shallower than the corresponding mean-field prediction, 
N(m) ~ mr'^l'^ . Note that the diffusion limited scaling 
exponent coincides with the reaction limited one when 
d = 2 which reflects the fact that c? = 2 is the criti- 
cal dimension for this model. In addition to having a 
different value for the scaling exponent of N{m), the sta- 
tistical properties of the diffusion limited regime are also 
structurally different: higher order correlation functions 
exhibit multi-scaling in the sense that they do not scale 
as powers of the corresponding single-point densities as 
they do in the mean field regime. 

In this article we focus in detail on this transition from 
mean field to fluctuation dominated regime in CCA. The 
fact that the parameter controlling the transition is the 
physical dimension of the model is inconvenient for both 
numerical exploration and physical applicability. We re- 



mark that the important ingredient for the application 
of mean field theory is not the physical dimension per 
se but rather whether or not the transport mechanism 
has the ability to break spatial correlations between par- 
ticles. We show that by replacing diffusive transport 
with Levy flights it is possible to continuously tune be- 
tween fluctuation-dominated behaviour and mean-field 
behaviour even in one dimension. The tuning param- 
eter is the characteristic exponent of the Levy distribu- 
tion underlying the jump size distribution for particle 
hops. The physical intuition is that correlations between 
particles are broken even in low dimensions if there are 
sufficiently many long range hops to effectively mix the 
clusters. We also explore the special role played in the 
forced case by the second order mass correlation function 
whose scaling exponent is strongly universal as a result 
of the so-called Constant Flux Relation (CFR) [10]. Its 
scaling exponent is, perhaps surprisingly, independent of 
both the spatial transport mechanism and the physical 
dimension and is always given by its mean field value. 

The outline of the article is as follows. We begin by 
defining the stochastic model used in studying CCA, as 
well as a brief overview of the statistical fleld theory used 
in deriving the Lee-Cardy equation for this model (Sec- 
tion 2). Then, we summarise the scaling properties of 
the multi-particle correlation functions or arbitrary order 
as applied to both reaction-limited and diffusion-limited 
regimes (Section 3). Using the Constant Flux Relation 
(CFR) we then derive the scaling invariance of the flux- 
carrying correlation function, demonstrating that this re- 
sult is independent of the dimension of the system and of 
the mechanism of spatial transport (Section 4). Imple- 
menting Levy flights for the one-dimensional model, we 
then characterise the suggested form of the mass density 
as a function of the Levy exponent, which we validate 
against numerical simulation together with taking direct 
measurements for the integrated forms of the mass cor- 
relation functions (Section 5). Here we verify the scaling 
independence of the flux-carrying correlation function as 
we alter the spatial transport mechanism. We finish the 
paper by determining a predicted form for the structure 
of the depletion zones as a function of the Levy expo- 
nent (Section 6), which we validate numerically before 
providing a conclusion. 



II. THE LEE-CARDY EQUATION AND 
STOCHASTIC SMOLUCHOWSKI EQUATION 

We now describe a simple stochastic model, first in- 
troduced without injection of particles in fllj and with 
injection of particles in [12|, which can be used to study 
the phenomenon of CCA. Discrete particles, each car- 
rying a mass index, m, occupy a d-dimensional square 
lattice, L, with no restrictions on the number of parti- 
cles occupying each lattice site. We take the mass index 
to be a non-negative integer multiple of a fundamental 
mass, mg. This means that all particles are composed of 
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clusters of a single basic monomer unit. The state of the 
system at time t is characterised by the set of occupa- 
tion numbers, N-^r^t, specifying the number of particles 
of each mass m at each point, x of the lattice. The occu- 
pation numbers evolve in time according to the following 
microscopic dynamical transitions which occur in contin- 
uous time with exponentially distributed waiting times: 

• Diffusion: A particle at any site may hop to an 
adjacent site. These hops occur at rate £), taken to 
be independent of particle mass. 

• Injection of monomers: A new monomer may ap- 
pear at any site on L. The injection of monomers 
takes place at rate J. 

• Aggregation: Two particles having masses mi and 
m2 which are at the same lattice point, x, may ag- 
gregate to create a single new particle having mass 
mi -I- rn2- This occurs at rate AA'x mi t (^x m2 1 ~ 
(^mi 7712 ) where A is the reaction rate, taken to be 
independent of mass, and the iVx m t factors encode 
the pairwise nature of the aggregation process (the 
Kronecker symbol, (5^^ , ensures that the rate is 
correct when two particles of the same mass aggre- 
gate). 

The "rates" in these rules refer to the corresponding wait- 
ing times. These rules define a Markov chain which can 
be easily simulated using the Gillespie algorithm [l^ . All 
numerical results presented in this paper have been ob- 
tained in this fashion. 

This model has a Master Equation whose solution can 
be used to compute averages with respect to the micro- 
scopic dynamics. Such averages, which we shall denote 
by are equivalent to averages with respect to the 
effective action of an associated statistical field theory. 
This effective action is obtained from the dynamical rules 
by the well known Doi-Peliti construction. For details 
of the procedure for this specific model see [M]. For 
many microscopic models describing interacting parti- 
cle systems, averages with respect to this effective action 
are further equivalent to averages with respect to solu- 
tions of a stochastic rate equation [HI known as the 
Lee-Cardy equation. The Lee-Cardy equation specific to 
our CCA model, first derived in [T7j, takes the form of 
a stochastic Smoluchowski equation. In the continuous 
limit it takes the form: 
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where ^xt is Gaussian white noise and A is the d- 
dimensional Laplacian. D, A and J are respectively the 
diffusion constant, reaction rate and mass injection rate 
obtained from an appropriate rescaling of the correspond- 
ing microscopic rates with the lattice spacing in the con- 
tinuous limit [14] . We shall denote averages with respect 
to the solutions of Eq. (fTl) by 

One of the most striking aspects of the Lee-Cardy 
equation is that the noise term is imaginary indicating 
that the field 0x m t cannot be literally interpreted as the 
density of particles, iVxmti however tempting such an 
identification may seem. Imaginary noise encodes the 
fact that the particles in the system are anti-correlated 
[TSl HHj. Physically, this anti-correlation comes about 
since the process of aggregation tends, on average, to 
produce a depletion zone around each particle. We shall 



discuss these depletion zones in detail in Sec. VI 



Careful study of the construction leading to Eq. ([T]) 
leads to the conclusion [14 that single point moments 
of 0x 771 1 correspond to factorial moments of the physical 
density, A'^xm*: 
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(2) 



Multi-point moments translate directly. For example, 

(0X mi t4^X 7712 t )i = {N^777i tA^xms t)L if mi ^ 7712. (3) 

These relations allow us to translate correlation functions 
of Eq. ([T]) which are most suitable for theoretical analy- 
sis into corresponding correlation functions in the origi- 
nal lattice model which are most suitable for numerical 
measurements. In relation to the interpretation of the 
field (/)xmt alluded to above, Eq. ^ tells us that 0xmt 
is equal to the density, N:x.77it, on average but differs in 
distribution. 

If the occupation numbers are large, we can, to leading 
order, ignore the distinction between moments and facto- 
rial moments on the right hand side of Eq. ([2| . Further- 
more, if we could neglect correlations between particles, 
higher order correlation functions could be factorised into 
products of densities: 
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This is the essence of the mean field approximation for 
CCA. If the density exhibits scaling with some physical 
quantity such as the mass or the source strength, it fol- 
lows from the mean field approximation that the higher 
order correlation functions exhibit simple scaling. That 
is to say that if the density has a scaling exponent, 71, 
then the corresponding scaling exponent of an order n 
correlation function is 7„ = 77 71 . 

We have already mentioned that the interactions be- 
tween particles resulting from aggregation provides a nat- 
ural mechanism for generating anti-correlations between 
particles. The extent to which the approximation Q is a 
good one therefore depends on the strength of these anti- 
correlations. It has been understood that in low dimen- 
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sions the large time, large mass behaviour is always dif- 
fusion limited resulting in anticorrelations between par- 
ticles which are sufficiently strong to invalidate the mean 
field approximation. Low dimensions in this case means 
d < 2 with the case d = 2 being marginal. One expects to 
observe a violation of mean field theory accompanied by 
the onset of multiscaling as one varies the dimension of 
the lattice L. This is indeed true with strong multiscal- 
ing observed in both the decaying [T7] and forced [TH HH] 
cases in d = 1 which mellows to logarithmic corrections to 
the mean-field scaling in d = 2. Triggering this transition 
from mean-field to non-mean-field behaviour by varying 
the spatial dimension is not very easy in practice. In 
Sec. |V] of this article we show that introducing long- 
range diffusion can break the anti-correlations between 
particles resulting in behaviour which is very reminiscent 
of mean-field behaviour even in one dimension. Before 
this however we will review the scaling properties of CCA 
with regular diffusion. 



III. SCALING PROPERTIES OF CCA WITH A 

SOURCE OF MONOMERS: REACTION 
LIMITED AND DIFFUSION LIMITED REGIMES 

A constant injection of monomers into the system al- 
lows for the creation of a non-equilibrium stationary 
state, which is characterised by a constant flux of mass 
from smaller to larger masses under aggregation. How- 
ever, since there is no upper mass limit defined on the 
system, the creation of arbitrarily large clusters can go 
on indefinitely. In this sense, the stationary state is only 
defined for some inertial range of masses. Any station- 
ary solution to the Stochastic Smoluchowski Equation 
describing the mass density must therefore be of Kol- 
mogorov type |S]. That is, these solutions exhibit scaling 
evocative of the Kolmogorov spectrum found in forced 
hydrodynamic turbulence, with a cascade defined by a 
constant flux of mass from smaller to larger clusters me- 
diated by the aggregation of particles. From this point 
on, we shall assume statistical homogeneity and drop the 
X index. 

The mass density has dimension [Nm] = L^^M^^. 
The dimensions of the physical parameters appearing in 
Eq. ^ are as follows [D] ^ L^T-^, [A] = L'^T-^ and 
[J] = M L~'^ . In order to describe a cascade, the 
stationary density should depend on the mass flux, J. 
There are two natural ways to construct an expression 
with the correct dimension using only one additional pa- 
rameter in addition to the flux: 
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The first corresponds to the reaction-limited regime in 
which diffusion plays no role. The second corresponds 



to the diffusion-limited regime in which the reaction rate 
plays no role since reactions are effectively instantaneous 
compared to the time required for particles to find each 
other by diffusion. 

If we average Eq. ([T]) in the statistically homogeneous 
case and make the mean field approximation, Eq. Q, 
then we obtain the standard Smoluchowski coagulation 
equation: 
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The exact stationary solution of this equation is O [7] 
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confirming that the reaction-limited scaling, Eq. ([s]), cor- 
responds to the mean field scaling. In one dimension, 
exact calculations [T^l showed that for the case of 
instantaneous reactions, A — >■ oo, 
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for m ^ 1, which agrees with the diffusion-limited scal- 
ing, Eq. ([6|, with d = 1. Furthermore the dynamical 
scaling properties of the time-relaxation to these steady 
states in both the mean-field and one-dimensional cases 
were computed exactly in [21]. It was subsequently 
shown in "22] that the diffusion- limited scaling, Eq. Q, is 
actually exact for all d < 2. The fact that the exponents 
in Eqs. ([5| and ([6| become equal when d = 2 reflects that 
the critical dimension for this system is 2. 

The presence of multiscaling in c? < 2 can be detected 
by measuring the higher order mass correlation functions: 



c„(mi, . . . ,m„,i) 
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The case rt = 1 is just the mass density as discussed 
above. In the statistically homogeneous case, these corre- 
lation functions do not depend on the spatial point, x, at 
which they are measured. For aggregation with source, 
c„(mi, . . . , m„, t) becomes independent of t as t — > oo 
and is a homogeneous function of mass for large masses: 

c„(/imi, . . . ,/im„) /i''"c„(mi, . . . ,m„). (11) 

In [H], the exponents 7„ were calculated perturbatively 
in e = 2 — d using dynamical renormalisation group: 
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and found to agree remarkably well with Monte Carlo 
simulations in d = 1 despite the fact that, formally, e = 1 
in this case. The epsilon expansion written in this way, 
with d appearing explicitly, demonstrates that the effect 
of fluctuations is to correct the simple scaling jn = nji 
(recall that 71 is given exactly by 71 = ^^p^ |22|), thereby 
demonstrating the presence of multi-scaling in the model. 
Full details are given in [T4. . Direct measurement of the 
exponents 7„ from numerics is quite difRcult for n > 1 
due to the strong fluctuations in the diffusion-limited 
regime. In practice, for n > 1, better statistics are ob- 
tained by measuring the integrated correlation functions, 

dm„_i c„(m, mi ... , m„_i). 

(13) 

Cn (to) can be expressed in terms of site occupation num- 
bers using Eqs. ^ and ([S]) and scales with exponent 
7„ — n -f- 1. Notice from Eq. (12) that 72 = 3, inde- 
pendent of the spatial dimension, to first order in the 
e-expansion. Furthermore, this is exactly the exponent 
one would expect from mean-field theory even though 



Cn{m) — / dmi 



Eq. ( 12 ) specifically describes the non-mean-field regime. 



In fact, 72 is exactly given by its mean- field value in all 
dimensions regardless of whether mean-field theory ap- 
plies or not. This is a specific case of a general property 
of many non-equilibrium stationary states which carry a 
flux of some conserved quantity: strong universality of 
the flux-carrying correlation function [10] . The fact that 
conservation laws determine the scaling of the flux-carry 
correlation function exactly is well know in turbulence 
as Kolniogorov's |-Law (see 23J for a review and wide- 
ranging discussions). We now study this important fea- 
ture in more detail. 



IV. CONSTANT FLUX RELATION 

We have already remarked that in the forced case, CCA 
reaches a quasi-stationary state characterised by a con- 
stant flux of mass from small to large masses. The Con- 
stant Flux Relation (CFR) [101 expresses the fact that 
this constant flux of mass exactly determines the scaling 
of the flux-carrying correlation function in the stationary 
state, a common feature of most driven non-equilibrium 
systems which reach stationary states dominated by a 
constant flux of a conserved quantity. 

The Lee-Cardy equation, Eq. ([T]), is equivalent to the 
underlying effective field theory, a fact which sometimes 
calls into question the usefulness of the Lee-Cardy for- 
malism. The derivation of the CFR, however, is enor- 
mously easier using the Lee-Cardy equation than it would 
be otherwise. At stationarity, and assuming spatial ho- 
mogeneity, Eq. ([1]) can be written 

= A / {(p{rn)(j){m — m'))^dm' 
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FIG. 1: Scaling universality of the second order integrated 
mass correlation function, C2(m) as given by Eq. (131, with 



mass, m, as shown for three different transport mechanisms. 
In comparison, scaling for the first and third order integrated 
mass correlation functions are shown in insets (A) and (B) 
respectively. 



where we have taken averages with respect to the noise 
term ^ and let m > ttiq. Let us assume that for large 
masses, the correlation function {(j){mi)(j){m2)) ^ takes the 
form 

C(TOi,m2) = (mima)"^^ f — ) , (15) 

\m2j 

where we have introduced the dimensionless scaling func- 
tion 'ijj. By the symmetry of C, ip has the property 
ip{x) = ip{x~^). This defines our flux-carrying corre- 
lation function. Re-writing Eq. (14) in the symmetric 
form 
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and applying the Zakharov transformations |24] 
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to the second and third integrals of Eq. ( 16 1 respectively. 



results in an expression in which each delta-function has 
the same argument. This gives the condition 
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which can only be satisfied for y := 2a; — 2 = 1, or equiv- 
alently, a; = |. Therefore, the flux-carrying correlation 
function must be of the form 



0(1711,7712) ^ (mirn2) ^ipi — ~l 

V' 



7712 J 



(19) 



for sufficiently large masses, which is a homogeneous 
function of degree -3. We have implicitly assumed that 



the corresponding integrals in Eq. ( 16 1 converge for this 
exponent. This can be demonstrated explicitly for the 
constant kernel case but, for general kernels must be ver- 
ified numerically |25| . An immediate corollary of this 
result applies to the specific case ra\ = ra2- using the 
factorial moments relationship in Eq. ([2]) we have that 
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Namely, that the stationary average density of pairs must 
scale with exponent -3 for sufficiently large m. 

In the mean field regime, higher order correlation func- 
tions factorise into products of densities as expressed in 
Eq. Q. It immediately follows from Eq. ^ that 72 = 3 
in mean field theory. Thus the scaling exponent, 72, is 
given by its mean field value in all dimensions. We note 
also that the dependence on the transport process dis- 
appears on averaging Eq. ([T]) so that the exponent 72 is 
strongly universal, as suggested by figure ([T]) . In the next 
section we consider CCA with Levy flights. The choice 
of the characteristic exponent, /?, of the Levy flights 
strongly affects the scaling of general correlation func- 
tions in the model. We shall see, however, that the value 
of 72 is insensitive to the choice of /?. 



V. ONE-DIMENSIONAL CCA WITH LEVY 
DIFFUSION 

CCA with Levy diffusion is of some interest in certain 
surface growth phenomena [571 HH] but our interest stems 
from the theoretical point of view that it allows us to 
mimic the transition to mean-field behaviour even in one 
dimension. From now on we focus entirely on the one 
dimensional case. 

We modify the model introduced earlier by allowing 
particles to make long range jumps. The waiting times 
between particle hops remain exponentially distributed 
with a fixed rate. Rather than simply hopping to nearest 
neighbours, however, the length of each jump is inde- 
pendently sampled from a probability distribution with 
a heavy tail. In our numerical experiments, we took this 
distribution to be a symmetric Levy distribution with 
scale factor 1 and scaling exponent /3 e (0,2]. Jump 
lengths were rounded to the nearest integer to remain 
consistent with our lattice formulation. Such hops are 
referred to as Levy flights [29] and should not be con- 
fused with the (often more physically realistic) case of 
Levy walks in which the waiting time distribution and 



FIG. 2: Comparison of the numerically measured exponent 
71 against the prediction derived from dimensional analysis. 



Eq. ( |23| l, as plotted against the Levy characteristic, /3 . Data 
was taken from stationary simulations on a lattice of size 10^ . 
The error bars on the exponents represent two standard devia- 
tions as estimated by bootstrapping a least squares estimator 
with the data measuring the integrated density, Ci(m), as 
defined by Eq. (1131), over 4 decades of m. 
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FIG. 3: Scaling exponents for the first three mass correlation 
functions as functions of Levy characteristic (3. The expo- 
nents 72 and 73 were estimated using a maximum likelihood 
estimator [26| for a power law distribution of pairs and triplets 
entering Eq. ( 13 1. Theoretically predicted values for both the 
mean field and diffusion limited cases are included for each 
exponent. Error bars correspond to two standard deviations 
calculated by bootstrapping the estimator. 



the jump length distribution are not independent to ac- 
count for the fact that for particles of finite mass long 
jumps should take more time. The probability distribu- 
tion, P{x, t) of the position of a particle starting at x = 
at i = and exhibiting one dimensional Levy flights with 
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exponent /3 satisfies a fractional diffusion equation: 
dP ^ d^P 

where the generaUsed (or anomalous) diffusion coefh- 
cient, Dp, has dimension L^T^^ and the fractional 
derivative is defined through its Fourier space represen- 
tation. See for details. For /3 = 2 the Levy dis- 
tribution reduces to the Gaussian distribution and we 
recover the case of regular diffusion. In this case, the 
mean square displacement of a particle grows linearly in 
time {x^it)) = Dt providing a natural way to define the 
diffusion length, Id = ^Dt. For /3 < 2, {x^{t)) is infinite 
for any finite time and we require a different approach 
to defining the anomalous diffusion length. One reason- 
able way to do this uses the fact that Eq. (21) has a 



self-similar solution describing the evolution of P{x, t) in 



which the self-similar variable is z 



\x\ 



[29]. It 



therefore makes sense, as one might expect from dimen- 
sional considerations, to define the anomalous diffusion 
length as Ip = {Dpt)^/^ despite the fact that the mean 
square displacement is divergent. 

One may repeat the dimensional argument outlined in 
Sec. |III| in the diffusion- limited regime taking the anoma- 
lous diffusion coefficient, Dp, instead of D and arrive for 
d = 1 at 



N„ 



oc 



J 



2 + 13 
' l + ;9 



(22) 



For /3 = 2 we recover the known diffusion limited scal- 
ing, Eq. Q. For /3 = 1, the scaling becomes that of 
the mean- field answer, Eq. (|8|, with Dp replacing the 
reaction rate, A. This suggests that at /3 = 1, the Levy 
flights have become sufficiently long range to break all lo- 
cal anti-correlations between particles. If /3 is decreased 
below 1, we do not expect that the scaling properties 
would change any further since particles are now effec- 
tively independent. Thus we have the following dimen- 
sional prediction for the scaling of the mass density in 
the presence of Levy flights: 



71 



2 

2+J3 



/? < 1 

1 < /3 < 2 



(23) 



Fig. |2] compares the results of numerical measurements 
of 7i with the predictions of Eq. ( 23 1 and indicates fair 
agreement. If we are correct in interpreting Eq. (23 1 as 



describing the breaking of correlations by the long-range 
hops then it should follow that multiscaling described by 
Eq. (12 1 should also be lost as /3 is tuned from 2 to 1 



and simple scaling should be restored for the higher or- 
der correlation functions. This is illustrated in Fig. [3] 
which shows the measured scaling exponents of the first 
three higher order correlation functions as a function of 
p. We see that within error, simple scaling is restored in 
one dimension as /? — >■ 1. This is not because we force 
the system to become reaction-limited by increasing the 



1.4 



1.2 



a 



m=1 

m=10 — B- 
m=100 — ^ 
m=1000 




50 



FIG. 4: The structure of the depletion zone as given by Eq. 
(271 in the case of regular diffusion. These plots have been 
normahsed against the single-point densities corresponding to 
each choice of mass, m. The inset shows plausible evidence 
in support of the predicted form given by Eq. (271 after an 
appropriate rescaling of the data by m}^^ . 



number of particles per site (the average number of parti- 
cle per site in all simulations shown here was about 0.8 so 
we remain in the low density regime). Rather the trans- 
port process breaks correlations between particles allow- 
ing the mean field behaviour to become observable again. 
It was pointed out to us that a somewhat similar transi- 
tion in the scaling exponents of CCA as a additional scal- 
ing parameter is varied was demonstrated within mean 
field theory in pf. There, an additional scaling parame- 
ter was introduced by the inclusion of exogeneous growth 
of clusters on top of the basic model dynamics. 

As an aside, we re mark that the dimensional argument 
leading to Eq. (23) suggests 71 = ^^p^ when d is left 
arbitrary. 



Setting 7i = 5 might lead one to conjecture 
that the critical dimension for CCA with Levy diffusion is 
dc = /3. While intriguing, this suggestion of a fractional 
critical dimension lacks the exact expression for 71 which 
allowed a convincing argument to be made for dc = 2 in 
the case of regular diffusion (/? = 2) [22] . 



VI. STRUCTURE OF THE DEPLETION ZONE 

The breaking of anti-correlations between particles by 
Levy flights should also be seen in the structure of the 
depletion zones surrounding particles. We now turn our 
attention to characterising this effect in one dimension. 
To detect the presence of depletion zones around particles 
in CCA we need to calculate the multipoint correlation 
functions between densities at different spatial points: 



Pn(xi, . 



t^xi nil i ' 



(24) 
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A closely related 2-point correlation function was calcu- 
lated exactly in one dimension for the case of infinite 
reaction rate in [52] and demonstrated clearly the ten- 
dency for clusters to anti-correlate in the diffusion-limited 
regime. It is natural to normalise such correlation func- 
tions with the mass density so that the ratio goes to 
1 when particles are uncorrelated. To avoid excessively 
clumsy notation, we suppress the mass dependence of the 
spatial correlations for the time being. 

If we forget about the mass dependence of the particles 
for the moment we are left with the coagulation process 
A+A — >■ A. One can study how the analogous correlation 
functions in this system behave. For the A + A ^ A 
system without source it was predicted in [3D| and later 
proven rigorously in |31| that in one dimension: 



5x 10'' 



4x 10" 



3x 10' 



^ 2x10 

d 

D 

1 X 10"^ 
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Regular Diffusion 
p=2.0 Levy diffusion 
p=1 .6 Levy diffusion 
p=1 .0 Levy diffusion 
Random flopping 




20 40 60 80 10( 



Pn(xi 




(25) 



where ~ denotes the short-distance asymptotic be- 
haviour. The first term on the right hand side is simply 
the nth power of the density. The second term accounts 
for correlations, producing multiscaling of the time-decay 
exponents for the A + A A problem in d = 1. The 
fact that p„(xi, . . . x„, t) vanishes as any two points are 
brought together reflects the fact that particles are anti- 
correlated. 

Since the mechanism for the generation of anticorrela- 
tions in the stationary CCA model is the same as for the 
A+A — A, one can use Eq. (25) to guess the correspond- 



ing formula for the CCA model in the stationary state 
[32] . In Eq. ( 25 ) the diffusion length for decaying coagula- 

The corresponding diffusion 



tion appears as Id = V Dt 
length for particles of mass m in stationary CCA is, from 
dimensional considerations, 



D m 



(26) 



Upon substitution into Eq. ( |25[ ) we arrive at the following 
suggestion for the asymptotic behaviour, now in mass, of 
the multi-point correlation functions in stationary CCA: 



Pn(xi 



. x„,m 



— TO 3 



n 



l<'i<j<n 



In the absence of a complete calculation, the plausibility 
of this argument should be assessed by comparison with 
numerics. This is done for n = 2 in Fig. |4] Eq. (271 im- 
plies that after rescaling with m^^^, p2(xi, X2, m) /pi{m)'^ 
for different values of to should collapse onto the same 
curve which vanishes linearly as the separation, |xi — X2I, 
tends to 0. Fig. |4] clearly shows that the depletion zone 
is larger for the heavier particles. The inset provides 
plausible evidence in support of the expected scaling. 

We may now ask what happens in the case of Levy 
flights. Fig. VI shows that as /3 is varied from 2 to 1, 



FIG. 5: The depletion zone as shown for different transport 
mechanisms without normalisation against the single-point 
densities, for a fixed choice of mass, m = 1000. As the Levy 
characteristic, P is varied from 2 down to 1, the anticorrelation 
void disappears, with the limiting cases of random lattice hops 
and regular diffusion also being shown. 



the depletion zone disappears. For comparison, the cor- 
responding correlation functions for the limiting cases of 
completely local hops and completely random hops are 
shown. The results for Levy flights lie somewhere in be- 
tween. Interestingly, the correlation between particles 
does not seem to vanish linearly with separation in the 
Levy case. See, for example, the main panel of Fig. |VI[ 
In the case of stationary CCA with Levy flights, the 
anomalous diffusion length for particles of mass to in one 
dimension is, dimensionally, 



Ids - 



Dp m 
J 



1+3 



(28) 



The fact that the correlation between particles vanishes 
linearly in the case of regular diffusion follows from the 
calculations performed in |30j and is connected to the 
fact that the spatial transport mechanism is diffusion. 



Dimensionally, any power of jx^ — x^ | //^j in Eq. (25) or 
Eq. (27 1 would be consistent so there is no reason to 



expect that this power should remain unity in the Levy 
case when the transport mechanism is changed: 



p„(xi, . . .x„,to) 



J\^ + I' _2±i3, 

— 1 TO ^ + f 



X n 

l<i<j<n 



^DpJ^^Y 1(1+15) 



(29) 



We put forward a scaling argument to fix the power 
of a. From the discussion of the Constant Flux Relation 
in Sec. 



IV 



we know that the second order correlation 
function at a single spatial point is theoretically propor- 
tional to the mass flux, J, in the stationary state regard- 
less of the transport mechanism. This fact is strongly 
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FIG. 6: Shape of the depletion zone as given by the predic- 
tion Eq. ( |30[ ) for the Levy characteristic /3 — 3/2. Each 
plot corresponds to a different choice of mass, m, and have 
been normalised against the corresponding single-point densi- 
ties. The inset provides plausible evidence for this prediction 
through the data collapse suggested by Eq. ( 30 1 . 



supported for the case of Levy flights by the numerical 
measurements presented in Fig. [3] If one accepts that 
this scaling remains true when the single point correla- 
tion function is split onto two nearby points, then one 
would expect that the power of J in Eq. ( 29 ) is unity 
for n ~ 2 (or equally, the power of m should be —3) 
which then requires that the exponent a should be equal 
to /3 — 1. This is consistent with Eq. ([27| when j3 — 2 



and consistent with the disappearance of the depletion 
zone when /3 = 1. For 1 < /3 < 2 we therefore propose 
the following shape for the depletion zone as measured 
by the two-point function: 



P2(X1,X2, 



TO 



4 + 2g 
1 + 3 . 



Xi - X2 



(30) 

This prediction is tested against the numerical measure- 
ments for /3 3/2 in Fig. VI The inset shows the mea- 

'Jpi{rnf' for several different 
i)/(i+/9) plotted as a func- 
According to Eq. (30), the data 



sured values of /Ci2(xi, X2 
values of to rescaled by wS^ 
tion of |xi — X2|''~"'^ 
should collapse to a single curve vanishing linearly for 
small separations. The collapsed curve seems to rule out 
the possibility of linear scaling as one might expect by eye 



from the main panel of Fig. VI and makes the suggested 
scahng, Eq. (30), seem plausible. Further investigations, 



preferably supported by analytic calculations, will be re- 
quired in order to draw a deflnitive conclusion. In partic- 
ular, it is important to properly characterise the "inter- 
nal" mass structure of correlation functions like Eq 
for different masses, particularly for n = 2. 
expected that the denominators of the \x. 
have mass dependences of the form {minijY^^ f{mi/mj 



(27) 



It might be 
— Xj I terms 



where f{x) is a homogeneous function of degree which 
has been ignored in our scaling arguments. The legiti- 
macy or otherwise of this neglect is likely to be somehow 
related to the question of the locality of the mass cascade 
(see [2S] for some discussion). This is a tricky issue in 
general and will require further work. 



VII. CONCLUSIONS 

To summarise, we have studied in detail, the scal- 
ing properties of one dimensional stationary CCA with 
Levy flights. We have demonstrated that the transition 
from mean field statistics to fluctuation dominated statis- 
tics usually observed as the physical dimension is tuned 
through two from above, can be mimicked in one dimen- 
sion by varying the characteristic exponent, /3, of the 
jump size distribution of the Levy flights. This is physi- 
cally reasonable since the introduction of long range hops 
provides a mechanism to weaken correlations between 
particles in the system and may erase them entirely if 
they are sufficiently frequent. Our predicted values for 
the scaling exponent of the mass density as a function of 
/3 based on dimensional arguments agreed well with nu- 
merical simulations of the underlying stochastic particle 
system. 

We also provided a direct demonstration of the strong 
universality of the scaling exponent of the mass-flux- 
carrying correlation function in the stationary state. It 
is independent of both the physical dimension and the 
value of P as expected from theoretical considerations. 

Finally we performed detailed investigations of the spa- 
tial structure of the depletion zones surrounding heavy 
particles in stationary CCA with both regular diffusion 
and Levy flights. Our results indicate that in the case 
of regular diflusion, the exact results of [30] describing 
the multiscaling of correlation functions in the decaying 
A + A ^ A model, have direct analogues for stationary 
CCA. In the case of Levy flights, our numerical studies 
suggest that the two-point correlation function measur- 
ing the anti-correlation between particles in the system 
vanishes with a non-trivial fractional power of the sepa- 
ration between the particles as this separation decreases 
to zero. Using our knowledge of the exact scaling for 
the flux-carrying correlation function we put forward a 
scaling argument that this power should be /3 — 1. This 
scaling is consistent with numerical observations but fur- 
ther efforts will be required to make it definitive. 

It is not immediately evident how to extend the meth- 
ods employed in [T31 [30] to account for Levy flights since 
the fractional Laplacian entering the propagator of effec- 
tive field theory seems difficult to obtain from the opera- 
tor representation of the master equation. A direct calcu- 
lation from the master equation seems equally daunting. 
Nevertheless we expect that some progress in this direc- 
tion could and should be made. This will greatly clarify 
the scaling arguments presented in this article. 

Several interesting extensions could be made to the 
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model in order to make it more physically realistic. One 
would be to introduce mass-dependent diffusion, for ex- 
ample, and allowing the diffusion rate to decrease as a 
power of the mass, D{m) ~ to"", as one would expect 
for real particles. Another would be to replace the Levy 
flights with Levy walks to account for the fact that phys- 
ical particles cannot instantaneously make long range 
jumps. We know that neither of these modifications 
would affect the scaling exponent, 72 of the two-point 
function which is fixed by the CFR alone. The value of 
7i and indeed 7„ for n > 2 are, however, expected to be 
sensitive to such modifications and would considerably 
complicate the picture. 

Another natural extension would be to ask similar 
questions about the structure of anticorrelations in the 
so-called charge model with Levy diffusion. The charge 
model [21] [33] is a closely related model in which masses 
(or charges) can take positive and negative values. When 
the rates of injection of positive and negative charges are 
equal, the flux J is zero on average so that the scal- 
ing arguments presented here fail. Instead of an exact 
conservation law for charge one instead has a statisti- 
cal conservation law for charge squared. This was used 
to derive the CFR (72 = 4) for the charge model in lOJ. 
Repeating the dimensional arguments leading to Eq. ( 22 ) 



in arbitrary dimension taking J to be the average flux of 
gives 71 — ^^^^ as a first guess for the charge model 

which gives the 



scaling. In d — 1 one has 71 



diffusion limited result 71 = 5/3 when /3 = 2 and the 
mean-field answer 71 = 2 when P = 1, both of which are 
known from |33| . If this dimensional prediction were sub- 
stantiated by an additional analytic or numerical study 
it could be used to generalise the remaining arguments. 
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